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ABSTRACT 

The completion of the GALEX All-Sky Survey in the ultraviolet allows activity measurements 
to be acquired for many more stars than is possible with the limited sensitivity of ROSAT or the 
limited sky coverage of Chandra, XMM, or spectroscopic surveys for line emission in the optical 
or ultraviolet. We have explored the use of GALEX photometry as an activity indicator, using as 
a calibration sample stars within 50 pc, representing the field, and in selected nearby associations, 
representing the youngest stages of stellar evolution. We present preliminary relations between 
UV flux and the optical activity indicator FH-g^ ^'^^ between UV flux and age. We demonstrate 
that far-UV (FUV, 1350-1780A) excess flux is roughly proportional to ^jjj^-We also detect a 
correlation between near- UV (NUV, 1780-2830A) flux and activity or age, but the effect is much 
more subtle, particularly for stars older than than ^ 0.5 — 1 Gyr. Both the FUV and NUV 
relations show large scatter, ~ 0.2 mag when predicting UV flux, ~ 0.18 dex when predicting 
"^HK' ^'^^ '^■^ '^^^ when predicting age. This scatter appears to be evenly split between 
observational errors in current state-of-the-art data and long-term activity variability in the 
sample stars. 

Subject headings: Galaxy: stellar content - stars: activity - stars: solar-type - ultraviolet: stars 



Introduction 



The highest regions of the solar atmosphere, above the visible light photosphere, are heated primarily by 
magnetic processes rather than by convection or radiative diffusion. The solar magnetic field is believed to 
be driven from the interface between t he radiative arid conv ective zones, although the details of this interface 
dynamo are not yet well understood (jCharbonneaul 120101 ). Since all Sun-like stars have a radiative and a 
convective zone, magnetic fields should be a generic feature of all Sun-like stars. As a proxy for magnetic 
fields associated with other stars, we can observe optical, ultraviolet, or X-ray emission from magnetically 
heated gas, as well as fiares, starspots, and other phenomena associated in the Sun with magnetic fields. 
These phenomena are collectively called stellar activity. 
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Stellar activity among Sun-like s tars decays on time scales of 10^ ^ yr (jSoderblom et al iPreibisch fc FeigelsonI 



2005t iMamaiek fc HillenbrandboOSi ) . much faster than the nuclear-burning time scale. This rapid eyolution is 
understood to be a consequence of magnetic braking slowing stellar rotation fe.g. lCharbonneau &: MacGregor 
1993 ) , in particular d ifferential rotation, and of differential rotation driving the generation of magnetic fields 
( Charbonneau 2010[ ). As a star ages it rotates more slowly, and thus its surface magnetic flux declines. 
A star with a lower magnetic flux cannot transport energy into its upper atmosphere as efficiently, and so 
the star's activity decays with time. Understanding the evolution of stellar activity with age is important 
because activity is one of the few age indicators that can be applied to main sequence stars. In addition, 
observing trends between stellar activity, age, rotation, mass, and metallicity can improve our understanding 
of the physics underlying activity in both other stars and the Sun. 

Because stellar activity is not a single phenomenon, it has many observational indicators. The most 
frequently used are the flux in the Ca II H and K lines at 3934, 68A, the soft X-ray flux, the Ha 6563A 
equivalent width, and the amplitude of a star's photometric variability. Of these, the Ca II lines are the 
most widely used because they have a number of advantages. They are resonance lines in the blue optical 
region, in which the non-thermal energy emitted by the outer stellar atmosphere can compete with the 
underlying thermal spectrum for detection. Other activity indicators accessible from the ground such as Ha 
or the Ca II IR triplet at 8498, 8542, and 8662A have much lower contrast, at least in solar-type stars, and 
are less efficiently observed. Stars from 0.5 to 1.0 in B — ^ (F8 V to K2 V spectral types) behave in similar 
fashion with respect to their Ca II emission. Stars as early as F5 V are also convective, but more weakly 
so, and stars later than K2 V show a different activity distribution in th e Ca II H and K lines, inc luding 
stronger emission than would be expected extrapolating from G dwarfs (jVaughan fc PrestonI 119801 ). The 
final advantage Ca II has over other activity indicators is that its measurements can be calibr a ted in terms 
of the Mount Wilson S-index, S*!^^, based on the original work of Vaughan fc Preston ( 1980l) . Noves et al 



(|1984l ) define the conversion of S'jyf^ into another standard index, -Rjjj^, representing the fraction of the 
star's bolometric luminosity emitted in the Ca II lines through chromospheric activity. These two papers 
together allow any Ca II flux measurement to be placed on a standard system, at least in principle. 

A star's uUraviolet emission lines, particularly the C II 1335A, C IV 1549A, and Mg II 2795, 02A 
features, have also been used as activity indicators. Like Ca II, they are high-contr ast resonanc e lines , 
with formation temperatures ranging from ^ lO"' K for Mg II to ^ 10^ K for C IV. ISimon et al.l (jl985h 
explored the relationship of UV emission to age among field stars, adopting ages based on lithium depletion 
calculations. They found that, between the ages of 1 Myr and 3 Gyr, UV lines tha t form at ^ 1 0'* K decay 
with an e-folding timescale of 2.7 Gyr, while ^ 10^ K lines decay faster, in 1.4 Gyr. iRibas et al.l (|2005l ) used 
a small sample of solar analogs to conclude that, between the ages of 100 Myr and 6.7 Gyr, the combined 
line and continuum flux in the 100-360A band decays as r"^'^" and that 920-1180A flux decays as r""-^^. 
So far, work in the UV band has been limited by the need to use space-based observatories, most of which 
have been designed for spectroscopy and can observe only small samples of stars. 



With the launch of the GALEX ultraviolet survey mission (jMartin et al.l 120051) and the release of the 
GALEX All-Sky Imaging Survey (AIS), it is possible to measure ultraviolet fluxes for many relatively bright 
stars (9 < F < 16 for G stars) entirely from archived observations. Such fluxes could be used to carry 
out an all-sky study of stellar activity of nearby stars, with more sensitivity than allo wed by ROSAT X-ra y 



observations, which detect only - 20% of nearby stars having 0.5 < B ~ V < 1.0 ICuillout et al.lll999l) 



Alternatively, GALEX fluxes provide a way t o select candidate ac tive stars for more detailed study. Such a 
selection has already been done for M dwarfs ( Shkolnik et ahlboil ). where activity produces an unambiguous 
UV excess. 
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In this paper, we explore the use of UV photometry as an activity indicator for Sun-like stars. In 
Section [2J we compare UV flux to the standard Ca II indicator, using a volume-limited sample of 

stars. In Section[3J we use young moving group and cluster members to look for a direct relationship between 
UV flux and age. In both sections we present empirical relations that predict ^jjj^ or age on the basis of 
UV flux, or vice versa. In Section |4l we verify that these relations do not suffer from a number of systematic 
errors. Finally, in Section [5] we discuss the implications of our results. 



2. UV Broadband Flux and R'j^j^ 

We explored the relationship between a star's flux in the GALEX FUV and NUV bands and its activity 
as measured by the -Rjjj^ chromospheric index using a uniform sample of spectroscopic data. In Section 12.11 
we present the ^jjj^ measurements, and in Section 12.21 the GALEX magnitudes. In Section 12.31 we present 
the method we used to fit relations between the UV magnitude and ^jjj^i examine the results of the fit, 
and construct UV-excess indicators, Rtj-it, defined analogously to Rtjtz- 



2.1. 



^HK ^^^^ 



A large-scale survey of HK emission in nearby G dwarfs was undertaken by one of us (D.R.S.) from 
1998 through 2002. The survey was done primarily using the Coude Feed telescope at Kitt Peak National 
Observatory to observe sources from the north pole to -40° declination. Camera 5 was used with grating 
RC250 in first order and with the F3KB CCD. A slit width of 250 gave good efficiency for the minimum 
resolving power of 4000 needed for the measurement. Stars below -40° declination were observed at Cerro 
Tololo Interamerican Observatory using the R-C spectrograph on th e 1.5 m telescope; the details of these 

CTIO observations included some 
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observations are the same as those reported in iHenrv et al. 
stars north of -40° to ensure cross-calibration. At both observatories various Mount Wilson standard stars 
spanning a broad range of color and activity levels were observed every run. Typical S/N levels were 
about 100 in the HK line cores but ranged over about 50 to 200 given changing observing conditions and 
approximate exposure times. The data were reduced using standard IRAF routines to take account of dark 
frames, bias frames, flat fields, and wavelength calibration exposures. 

The -Rjjj^ mea surements are on the Mount Wilson system, defined by the spectrograph built by 
Vaughan fc Preston ( 1980l ) with triangular band-passes centered on each of the H and K lines of width 
of I.09A and two 20A continuum bands on either side of the HK region. These filters are used to measure 
•^MW' ^hich is the ratio of the fiux in the cores of the H and K lines to that in the continuum. A correction 
for the star's spectral type, using the star's B — V color, is then made to yield ^jjK' ratio of flux in the 
cores of the H and K lines relative to the star's bolometric fiux. A c orrection for photo spheric flux in the 
HK bandpasses, also based on i? — V^, is then made to produce Rjjj^ (|Noves et al.lll984() . 

Several stars in the sample were observed multiple times; the average scatter in logi?jjj^ for these 
stars is 0.0813 dex. Our estimate of t he scatter is consistent with previously publis hed estimates, which are 
typically 0.05 dex in the i^^dex (jWright et al.ll2004Hlsaacson fc FischenbOlOl ). The index shows 

less scatter than the ^jjj^ index because its flux includes a constant photosphere contribution, but we can 
correct for this distinction by multiplying the 0.05 dex figure by the median ratio of the Ca II index with 
and without the photosphere, ^^H K follows that previous studies have seen scatter in -Rjjj^ 

of ^ 0.085 dex, similar to our own figure. ISoderblomI (|1985[ ) suggests the systematic measurement errors in 
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^HK ^""^ ^ '^■^^ ^ 0.07 dex of the scatter from intrinsic variability. 

The goal was to observe a sample consisting of all G dwarfs within 50 pc. A star was considered to be 
a G dwarf if it fell within 1 magnitude of a solar-composition zero-age main sequence (ZAMS) according to 
astrometry and photometry from the Hipparcos Catalog. Specifically, the star had to have a B — V color 
from 0.50 to 1.00, inclusive, a Hipparcos parallax of 20 mas or larger with no regard paid to parallax error, 
and an absolute V magnitude (determined from the V magnitude in the Hipparcos catalog and the Hipparcos 
parallax) below a line constructed to be 1.0 magnitude above a ZAMS. A close approximation to this is a 
straight line drawn from My = 3 a,t B — V = 0.5 to My = 5.5 a,t B ~ V — 1.0. As one would expect, 
the sample is distributed isotropically across the sky, except for one region containing excess stellar surface 
density from the Hyades cluster. 

The sample observed is not a pure one in several regards. Ideally one would like to be able to define and 
then observe a complete, volume-limited sample of stars. In reality several effects are taking place. First, 
the Hipparcos sample is imperfect. The Hipparcos Input Catalog attempted to be magnitude- limited to a 
V magnitude near 9 (essentially the HD catalog). However, for technical reasons the mission could observe 
a maximum number of stars in any one region and that means that some fairly bright stars are missing and 
that the magnitude limit has a ragged and inconsistent edge. Neither of these effects is likely to bias our 
sample because they were applied for strictly technical reasons with little regard for source type. 

Second, the overall magnitude limit of the HIC does impose a bias, both in age and multiplicity. The 
Sun would have y 9 if it were 50 pc away. Solar mass and composition stars younger than the Sun are 
not included beyond 50 pc due to their lower luminosity. Lower-mass stars suffer a larger age bias with 
a distance horizon for completeness that gets lower fast for cooler stars. However, we have tested for and 
found no detectable bias in the mean parallax with stellar color (see Section 1431) . suggesting that this is not 
a major effect. Binary systems can be brighter than single stars and so are more likely to be included in the 
HIC. 

Third, the sample is only as pure as the data used to define it. This is generally not a problem because 
the stars in the sample are bright and well-studied, but sometimes the colors or magnitudes may be wrong. 
Also, the Hipparcos parallaxes have errors, and so there are some stars in the sample that are further away 
than we think and perhaps a few with underestimated parallaxes. 

Fourth, we found obvious non-dwarfs during the course of the observing program. As noted, we paid 
no regard to parallax error in constructing our sample. One consequence is that the sample contains more 
distant stars that are non-dwarfs and that have large parallax errors, appearing to place them within our 
domain. Some of these appear to be cool white dwarfs or reddened A stars. We removed these contaminants 
from the sample, leaving 2680 stars. 

2.2. GALEX Data 

GALEX has a lower spatial resolution (~ 5") than Hipparcos, allowing cases where flux from a single 
GALEX source must be arbitrarily divided between two optical sources. In addition, double stars unresolved 
by either GALEX or Hipparcos could show discrepant colors that would confuse the analysis. To avoid such 
systematics, we removed known binary stars that were too close for GALEX to resolve, or whose separation 
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we could not determine, by dropping stars meeting the following criteria in the Hipparcos catalog: 

• Either the proximity flag is set, there is more than one stellar component in the Hipparcos entry, or 
the multiple systems flag is equal to C, G, or O, and 

• The component separation is less than 10" or not specified 

Of the 2680 stars in the original sample, we removed 19% as potentially unresolved binaries. 

We searched for counterparts from the GALEX All-Sky Imaging Survey (AIS) GR5 data release^ for 
each star remaining in the -Rjjj^ sample. Many of the stars lay outside the survey's sky coverage, either 
because they were bright enough to damage the GALEX detectors or because they were located too close to 
the Galactic plane, where the number of stars poses a similar danger. We required that GALEX sources be 
within 7" of the Hipparcos position, corrected for proper motion. The large search radius, compared to the 
GALEX astrometric accuracy of 0.5", was needed because the GALEX PSF is considerably degraded for 
bright stars. We also required that sources be within 0.5° of the field center to minimize uncertainties in the 
PSF and the detector response. We did not filter the GALEX sources by artifact flags, as the flags (such as 
"proximity to a bright source" ) are all either invalid or irrelevant for bright stars. After removing stars that 
were not covered by the AIS, we were left with 1204 stars, of which all but one were detected by GALEX. 
HIP 48141 was not detected because it was observed only in the FUV band, where we are not sensitive to 
cooler stars. 

We used aperture photometry provided as part of the GALEX GR5 source catalojl, using a 17.25" 



radius . We added an 0.07 mag aperture correction in both the FUV and NUV, following iMorrissev et al 



(j2007l ). The median photometric error in the matched sample is 0.006 mag in the NUV, but only 0.13 mag 



in the FUV because many of the sources are faint in the FUV. We show the matched sample in Figure [T] 

Many of the stars observed by GALEX were significantly saturated in the NUV. Since we cannot yet 
reliably correct the fiuxes of heavily saturated GALEX sources, we could not use these stars in this paper. 
Instead of making a cut directly on NUV, which would bias the sample towards intrinsically UV-faint stars, 
we required that V + 6.4:6{B — V) > 12.8. This cut, whose slope we derived from a linear fit to the NUV — V 
and B — V colors, ensured that most stars in our sample had NUV > 13.9, the magnitude below which 



saturation lowers the apparent flux by 10% or more (jMorrissev et al.ll2007l) . but at the cost of removing most 



stars with B — V < 0.7 from our NUV sample. We illustrate the efficiency of the cut in Figure [TJ 

Many of the stars were not detected in the FUV. We constructed formal detection limits for these stars 
by finding the smallest flux that would give a 3cr flux measurement. We used an SNR cutoff, rather than 
trying to estimate the threshold below which the GALEX pipeline no longer identifles sources, to avoid any 
stochastic behavior in the source extraction algorithm at faint levels and to avoid including S/N ^ 2 sources 
whose formal flux measurement provides almost no information about the source flux. We modeled the 
photometric noise as \/ AB + F, where A is the aperture area, B is the sky background (obtained, via the 
source catalog, from the GALEX pipeline's smoothed background maps), and F is the source flux. According 



^The official Hipparcos designations of the fields we used are as follows; proximity flag is H2, number of components is If58, 
multiple systems flag is H59, and component separation is H64. 

• 'jhttp: //galex. sts ci . edu/GR4/ [sic] 

• ^http: //galex. stsci . edu/casjobs/ 1 
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to this model, to be detected at signal-to-noise ratio S a star must have a flux 

F=1{S^ + VS'^ + 4S^AB^ (1) 

If a star was not detected, we quoted its flux as an upper limit at the value given by Equation [1] with S* = 3. 
If a star was identified by the GALEX pipeline, but with a recorded flux less than that given by Equation [I] 
we also reported it as an upper limit. This ensured that our detection limits were consistent across the entire 
sample, whether or not SExtractor was able to identify a source at lower SNR. 



2.3. Results 

2.3.1. Color- Color Locus 

To illustrate the effect activity has on stellar SEDs, we plot in Figure [2] UV — vs. B — V plots, 
with stars segregated by -Rjjj^, and in Figure [3] J7T^ — V vs. -Rjjj^ plots, with stars segregated hy B — V. 
Because all stars in the sample are within 50 pc, we assume no reddening. A clear trend appears in the 
stars' FUV — V color, in that more active stars with higher values of -Rjjj^ consistently have bluer FUV-V 
colors. The NUV — V color shows a similar, but much weaker trend; only for stars with B — V > 0.9 can 
we clearly distinguish between active and inactive stars. This behavior is not surprising. The photosphere 
contributes a larger fraction of the flux in the NUV than in the FUV over our investigated color range, 
0.5 < B — V < 1.0, as we expect from T^/f < 6200 K photospheres and as we confirm in Section [2.3.3l Since 
any variation in stellar activity is diluted by the photosphere flux, we expect to see a more subtle effect in 
the NUV. 



2.3.2. Nonparametric Fits 



In this section we investigate whether UV and optical activity indicators trace one another. Our goal is to 
convert between UV flux and -Rjjj^, and use one as a predictor for the other. In the absence of a quantitative 
physical model for activity indicators in either the UV or the optical, we fit polynomial expansions to the 
data. To pr edict UV acti vity from Ca II activity, we used an expectation-maximization (EM) algorithm as 
outlined by lAitkinI (|l98lh to incorporate upper limits to the FUV flux. Predicting ^jjj^ from UV flux was 
more challenging. Since we know of no algorithm for handling limits in independent variables, we were forced 
to restrict our FUV sample to avoid undetected stars. In analogy to our procedure for avoiding saturated 
NUV stars, we removed all stars from the FUV sample with V + 12.03(5 — V)> 15.5 when fitting for -Rjjj^ 
as predicted by UV flux. As a result, while we predict FUV flux for all stars with 0.5 < B — V < 0.9, our 
prediction of ^jjj^ from FUV flux holds only if i? — V" < 0.7. 

For both the predictor of UV flux and the predictor of ^jjj^i we chose the number of terms in the 
expansion that minimized the leave-one-out cross-validation score CV — (l/?i) "^jjiui — f(-i){xi)Y ^ where 
D is the set of stars detected in the UV, / is our estimate of the true relation f{X) = EiY\X), and is 
the fit we would have obtained leaving out the data point {xi,yi). Using CV as a goodness-of-fit statistic. 



and minimizing it, approximately minimizes the expected mean square error E 



independent of the functional form of the fit or the error distribution around it (jWasserman 



m^e.g.). 



We recognize that the procedure of lAitkinI (jl98l[ ) has a number of limitations. It cannot be extended to 
known measurement errors, and in fact implicitly assumes that the scatter is intrinsic to the system being 
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studied. Unfortunately, to our knowledge there is no good algorithm that both avoids similar assumptions 
and applies to models more complex than a straight line. In light of such concerns, we test the robustness 
of our fits in Section |31 

We found the following fits between a star's UV — V color and its R'^iK ^^1^®- -^^^ ^'-'^'^ relations require 
a. B — V color to control for the dependence of UV — V color on spectral type. After each equation we give 
the region of UV — V vs. B — V oi -Rjjj^ vs. B — V parameter space where it is valid. 

FUV~V = 12.30- 3.95(logi?jjj^ +4.5) - 2.94(logi?^j^ +4.5)2 

+ 2.22(5 -V- 0.8) - 9.11(5 -V- 0.8)(log + 4.5) 

- 12.0(5 -F- 0.8) (log i?^j^+ 4.5)2 

- 17.6(5 -V- 0.8)2 _ 3 _y_ o.8)2(log R'^^ + 4.5) 

- 24.3(5 -l/-0.8)2(logi?jjj^+ 4.5)2 (2) 
(0.5 <B -V <0.9 and 1.7(5 - V) < log + 6.62 and - 5.2 < log < -4.2) 



log5jjj^ = -4.25 - 0.44(i^[/F - V" - 12) 

+ 3.48(5 -V- 0.8) - 0.35(5 -V - 0.8){FUV - - 12) (3) 
(0.5 < 5 - F < 0.7 and 3.4 < FUV - F - 12.0(5 - V) < 5.2) 



NUV-V = 6.19- 0.87(logi?^j^ +4.5) - 0.93(logi?^j^ +4.5)2 

+ 6.54(5 -V- 0.8) - 2.73(5 ~V~ 0.8)(log5jjj^ + 4.5) 

- 2.89(5 - F-0.8) (log +4.5)2 (4) 
(0.7 < 5 - F < 1.0 and - 5.2 < logi?W < -4.2) 



log5jjj^ = -5.17-0.29(W\/-T/-7) + 0.64(iVJ7F- V^-7)2 
+ 0.28{NUV -V-7)^ + 4.06(5 -V - 0.8) 

- 2.68(5 -V- 0.8){NUV -V -7) 

- 2.15(5 -V- 0.8){NUV -V -if 

-Q.bl{B -V ~{).8){NUV -V -if (5) 
(0.7 < 5 - 1/ < 1.0 and 0.7 < NUV -V - 6.5(5 -V)< 1.6) 

Like all physically unmotivated fits. Equations [2][5] are best thought of as interpolations over the region of 
parameter space where we have data, rather than as functions whose form has physical significance. 

The RMS residuals around Equations [2]l5] are 0.33 mag, 0.15 dex, 0.21 mag, and 0.18 dex, respectively. 
Propagating the formal photometric errors and the observed scatter in log-Rjjj^, the expected residuals 
around Equations [2][5] are 0.28 mag, 0.10 dex, 0.16 mag, and 0.28 dex, respectively. Since the error propaga- 
tion makes a variety of assumptions, such as Gaussian errors, which are only approximately true, we do not 
expect a perfect match between the expected and actual residuals. 

The expected residuals have a large contribution from the photosphere color B — V, because of the steep 
dependence on UV photosphere flux on effective temperature, and a comparable contribution from scatter 
in ^jjj^- For Equations [2] and [3l scatter in ^jjj^ contributes ^ 60% of the error budget, with uncertainties 
in FUV — V and B — V contributing ^ 30% and ^ 10%, respectively. The expected residuals around 
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Equation |4] are dominated by uncertainties in B — V (84%), plus a 14% contribution from scatter in R'^j^- 
The residuals around Equation [S] are evenly split between uncertainties in B — V and ^jjj^ (45% and 54%, 
respectively). If ^ 75% of the variance in the ^jjj^ measurements is from intrinsic variability fSection l2.ip . 
then variability in R'^j^ accounts for 40-45% of the residuals around Equations [21 [21 and[5l 



2.3.3. Normalized Excess Fluxes 

Equations [2][5] are purely empirical results, unnormalized by the stellar photosphere and independent of 
any theoretical expectations of how much UV flux an active star will produce. However, it is often useful 
to divide the UV flux into a fixed photosphere component and an excess associated with stellar activity. 
Therefore, we also used Kurucz photosphere models to calculate broadband UV excess activity indices 
Rjjyj in analogy to the spectroscopic activity index ^jjj^, as follows. 



The GALEX magnitudes are defined on the AB system (|Okdll974r ). allowing us to estimate the total 
flux observed in the GALEX bands as 

/yy = (1 erg s-i cm-2 Hz"!) x io-"-4(c/v+48.60) ^ Aiy-^y (6) 

where Au is the frequency range corresponding to the wavelength range 1350-1780A (for FUV) or 1780- 
2830A (f or NUV). The adopted bandpasses represent the wavelengths where the effective area falls to 10% of 
its peak ( Morrissev et al.ll2007 ). rounded to lOA precision for notational convenience. While the bandpasses 
are somewhat arbitrary, we show in Section 14.41 that they make the fluxes we calculate insensitive to the 
color of the source at the ^ 10% level. 

We found bolometric corre ctions and photosphe r e UV fluxes to each star in our sample, using the 
spectral type fitting method of Kraus fc Hillenbrandl (2007). We constrained the fits using all available 
photometry, namely the Hipparcos B and V magnitudes and J, H, and K magnitudes from the Two-Micron 
All Sky Survey. Table [T] summarizes the run of empirical magnitudes in these five bands and bolometric 
magnitudes along the main sequence, based on lKraus fc Hillenbrandl (|2007l) Table 5. To fit photosphere UV 
fluxes we also needed UV magnitudes along the main sequence, which we obtained from solar-metallicity 
Kurucz "ODFNEW" modelf0 for each spectral type. We matched the Kurucz grid to each spectral type by 
using the masses, effective temperatures and bolometric magnitudes in lKraus fc Hillenbrandl (|2007l) Table 5 
to infer values of iogg at each spectral type, then linearly interpolated Kurucz models to the (Tg-g-, log g) 
pair. Since UV and optical magnitudes are normalized on two different systems , we us ed radii inferred from 
the bolometric magnitudes and effective temperatures in Kraus fc Hillenbrandl (2007) Table 5 to place the 
model fluxes on an absolute scale. 

We then calculated absolute FUV and NUV magnitudes for each interpolated model by integrating 
the photon flux over GALEX effective area curves provided by A. Gil de Paz, and dividing by photon fluxes 
i ntegrated over a GAL SPECH spectrum (version mod_002) of LDS 749B, the primary GALEX calibrator 
(jMorrissev et al.l 120071 ). Working with photon fluxes allowed our synthetic magnitudes to automatically 
incorporate color corrections, which can grow up to ^ 0.2 mag for K stars. We also calculated absolute B 
and V magnitudes using response curves from iBesselll (|l990bf ). a CALSPEC s pectrum of Vega (stis_005), 
and the Landolt magnitudes oi B = 0.02, V = 0.03 for Vega (jBesselll Il990al ). This ensured that, while 



' jhttp : //kurucz ■ harvard ■ edu/j 

' ^http : //www, stsci . edu/hst/observatory/cdbs/calspec .htrnT] 
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these B and V magnitudes could not be considered color-corrected in the same sense as the FUV and NUV 
magnitudes, all four magnitudes were calculated consistently. 

When the FUV and NUV magnitudes calculated from each model were matched to the spectral type 
the model was intended to represent, they produced a curve that did not match our observed stellar locus 
in color-color plots, as shown in Figure [5] In particular, the predicted photosphere NUV magnitude was 
typically larger than the observed (photosphere plus activity) NUV magnitudes. Investigating, we found that 
our synthetic ^ vs. B — V relation closely matched the data, while the B—V vs. T^g- did not, suggesting that 
the problem lay in the effective temperature adopted for each spectral type. We created the UV magnitudes 
we list for each spectral type in Table [T] by linearly interpolating the model UV — V vs. B — V relation to 
the B — V observed for each spectral type, and adding the observed V magnitude. Figure [5] shows that these 
corrected magnitudes appear much more consistent with our observations, particularly in the NUV. 

Having constructed Table [1] we estimated bolometric , phot osphere FUV, and photosphere NUV mag- 
nitudes for each star as described in iKraus fc Hillenbrandl |2007t ). We then calculated the activity indices 



r rphot r „phot 

_ JFuv - /FUV _ mvv - /nuv .7^ 

analogous to ^jjj^- Here, /uv ^^"^ inferred from the observed FUV or NUV magnitude, as defined at 

the beginning of this section, f^jy^ is the flux inferred in the same way from the estimated photosphere FUV 
or NUV magnitude, and /j^qJ is the estimated bolometric flux. We propagated errors on Rjjy neglecting 

any correlation between f^^^ and f\,Q\', even if the two estimates were perfectly correlated, the correlation 
term would constitute only ^ 2% of the error budget. The median uncertainty in logi?puy was 0.082 dex, 
80% of which was propagated from uncertainties in /puV" "^^^ uncertainty in log-Rj^uy was 0.100 dex, 

^ 95% from uncertainties in /^^jy- 

We compare the R\jy R'^^ indices in Figure ID ^?j.^uv ^^^"^^ correlation with R'^^ , and a 
mean of —2.0 x 10~^ ± 3.1 x 10~^, consistent with zero. Since we do see a correlation between NUV — V 
with -Rjjj^ , we believe our lack of a correlation for ^jyfjjy represents the uncertainties in the photospheric 

contribution /pj^y overwhelming any variation in NUV flux caused by activity, ^ puy is well correlated 



with ^jjj^- We fit a line to the logs of both activity indicators using the procedure of lAitkinI (jl98l[) . getting 
the relation 

log i?puy = (0.98 ± 0.05) log -I- (-0.53 ±0.25) (8) 
We compare our estimates of ^pyy to previous work on UV stellar activity in Section 



3. UV Broadband Flux and Stellar Age 

The volume-limited Hipparcos sample of Section [5] is uniform and assumed to be extinction- free, but 
because it consists of field stars we do not have precise age estimates for individual members, nor do we 
probe ages much below 625 Myr, the age of the Hyades. To extend the study to younger ages, and to 
explore the dependence of UV flux on age itself rather than on an empirical proxy for age, we constructed 
a second sample consisting of members of nearby clusters and moving groups. We present this sample and 
the matched GALEX data in Section IXTl In Section [5?^ we present our relations between UV magnitude 
and age. Since our procedures for the moving group sample parallel those we used for the R'^^ sample, we 
summarize them here and refer the reader to Section [5] for details. 
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3.1. Data 



We selected members of the TW Hydrae (8 Myr old), (5 Picto ris (12 Myr old), Tucana/ Horologmm 
(30 Myr old), and AB Doradus (50 Myr old) movi ng groups from Zuckerman &: Song (|2004l ). as well as 
mem bers of the clusters Bla nco 1 (120 Myr old) from iMermilliod et alj (j2008( ) and the Hyades (625 Myr old) 
from Ferryman et al. ( 19981 ). All ages are from the same papers as the membership lists. We could use only 
two clusters as age benchmarks because the vast majority of nearby clusters, such as the Pleiades, IC 2602, 
Alpha Persei, and Coma Berenices, have members bright enough to trip the GALEX safety limits and so 
could not be observed. More distant clusters, such as M 67 and NGC 188, are too distant for GALEX to 
detect members cooler than late G type, so these clusters would produce samples that have little overlap 
in spectral type with the nearer groups. We did not restrict membership to stars with 0.5 < _B — < 1.0, 
as in Section [21 so compared to the Hipparcos sample the moving group and cluster sample has many more 
low-mass stars and a handful of higher-mass ones. 

As in the Hipparcos sample, we removed known binary stars that were too close for GALEX to resolve, 
or whose separation we could not determine, to avoid systematics associated with blended sources. Since 
the data set was inhomogeneous, the criteria for removing binary stars varied by group: 



Hyades: iPerrvman et al.l (|l998l ) considered only Hipparcos stars, so we dropped stars following the same 
rule as in Section! 



• Either the proximity flag is set, there is more than one stellar component in the Hipparcos entry, 
or the multiple systems flag is equal to C, G, or O, and 

• The component separation is less than 10" or not specified 

Blanco 1: we dropped stars flagged by Mermilliod et al. ( 20081) as resolved double stars or spectroscopic 
binaries 



Moving groups: we used the SIMBAD database to look up literature on each star, and dropped stars with 
companions within 10". 



Of the 425 stars in the original sample, we removed 26% as possible unresolved binaries in GALEX obser- 
vations. 

Many moving group members, particularly the lower-mass stars, do not have high quality optical pho- 
tometry. Instead, we used infrared photometry from the Two-Micron All Sky Survey (2MASS) for our 
analysis, even though using infrared rather than optical colors increases the scatter around our relations (see 
Appendix Each moving group or cluster member had exactly one 2MASS counterpart within 3", so the 
matching was straightforward. The median errors on the J, H, and K magnitudes in the matched sample 
were 0.02, 0.03, and 0.02 mag, respectively. 

We searched for counterparts within 7" of each star in the GALEX All-Sky Imaging Survey (AIS) GR5 
data release. As with the Hipparcos sample, many of the stars could not be observed by GALEX because of 
the observatory's brightness limits. GALEX coverage proved a much more serious restriction than for the 
Hipparcos sample. Only 99 of the 313 single stars in the sample were observed by GALEX. The most badly 
affected group was Blanco 1, where only 6 outlying stars were observed out of 49 single members; the more 
centrally located members were too close to bright stars. Of the 99 targets three (TWA 26, TWA 28, and 
HD 89744 B), all brown dwarfs, were undetected by GALEX. 
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We processed GALEX sources in the same way as in Section [2^ we considered only sources within 0.5° 
of the field center and used corrected aperture photometry from SExtractor. We avoided NUV-saturated 
stars by only considering stars with J + 7.85( J — K) > 10.3 for NUV fits, and we constructed upper limits 
for missing fluxes from the exposure time and the local sky background using Equation [TJ The median FUV 
error for detected sources was 0.14 mag, while the median NUV error was 0.01 mag. 



3.2. Results 

We plot in Figure[7][/F— J vs. J — K plots, with stars segregated by group or cluster membership, and 
in Figures [BJH] C^t^ — J vs. age plots, with stars segregated hy J — K . There is a clear correlation between 
NUV — J and age, in the sense that younger stars have bluer NUV ~ J colors at fixed J ~ K. The FUV 
data, on the other hand, are confused by the large number of non-detections. We believe we see a stronger 
trend between NUV flux and age here than we saw between NUV flux and -Rjjj^ in Section 12.31 because we 
are considering younger stars, whose activity levels drop more rapidly than in the volume-limited sample. 

Although these stars are generally more distant than those presented in Section [2l exti nction is still 



negli gible. The most heavily reddened group in our sample, Blanco 1, has E(B — V) = 0.01 6 (jCargile et al 



20091 ). which implies E{FUV - J) ~ E{NUV - J) - 0.1 based on the Cardelli et al.l (|l989r ) reddening law. 
Since both the trends and the scatter in Figures [7] and [8] are much larger than 0.1 mag, we ignore foreground 
extinction. 



We fit polynomial curves to the NUV data following the same procedure as in Section [2731 We also tried 
to fit the FUV data, but the results were quite poor - because many of the stars in this sample were cooler or 
more distant than those in the Hipparcos sample, we did not have enough FUV detections to meaningfully 
constrain the relationship. 

We found the following fits between a star's NUV — J color and its age. All four relations require a 
J — K color to control for the dependence of NUV — J color on spectral type. After each equation we give 
the region of NUV — J vs. J ~ K oi log age vs. J — K parameter space where it is valid. 

NUV - J = 8.77 + 0.79(log(age/yr)-8.0) + 8.91(J- A'-0.6) 

+ 0.87( J -K- 0.6)(log (age/yr) - 8.0) - 3.57(J -K- 0.6)^ 

- 3.92( J -K- 0.6)2(log (age/yr) - 8.0) (9) 
(0.4 < J - if < 0.9 and 6.9 < log (age/yr) < 8.8) 

log (age/yr) = 7.53 -I- 0.80(iV[/F - J - 8) 

- 7.69( J -K- 0.6) + 0.12(J -K- Q.Q){NUV - J - 8) (10) 
(0.4 < J - /s: < 0.9 and 2.7 < NUV ~ J ~ 7M{J - K) < 5.3) 

Equations l9l and fTOl have residuals of 0.46 mag and 0.39 dex, respectively. The expected residuals around 
these two fits are 0.63 mag and 0.43 dex from propagated errors in NUV — J and J — K. We did not consider 
scatter from UV variability, although the NUV magnitudes of M stars at the ages of the Pleiades and Hyades 



vary by ~ 1 mag (jBrowne et al.ll2009l ). and one might expect variability to be detectable against the brighter 
photospheres of G and K stars as well. Errors in J — K dominate the expected residuals, accounting for 97% 
of the variance. We do not consider the effect of age errors, since they are highly correlated (i.e. our age 



- 12 - 



estimates for all members of a group are in error by the same amount) and to first order should affect only 
the fits and not the scatter around the fits. As with the Hipparcos sample, we can account for the size of 
the residuals from propagated errors alone, particularly from poor J and K photometry for stars saturated 
in 2MASS. 

We also used the procedure of 12.3. 31 to estimate bolometric, photosphere FUV, and photosphere NUV 
magnitudes for our moving group members using their 2MASS photometry but no i? or data. The estimates 
were possible only for the bluer stars (inferred spectral type K5 or hotter) , as the limited temperature coverage 
of the Kurucz grids did not allow us to find photosphere UV fluxes for cooler stars. We found a median error 
in log-Rpyy of 0.089 dex, split evenly between propagated uncertainties from /puv ^'^"^ from /pjjy- The 
median error in log/Jj^^uy was 0.26 dex, 95% of which came from uncertainties in /^yy- 

We plot R'jjY vs. age in Figure (TU] As in Figure IHl -Rpyy shows a trend with age, albeit with several 
orders of magnitude of scatter. There is a hint of a trend in ^Js^uy sls well, but as in Section 12.3.31 the 
calculation of -Rj^yy is dominated by noise. We were able to fit the evolution of the FUV fiux as 

log i?^UV = ("0.42 ±0.18) log (age/yr) + (-1.27 ±1.50) (11) 

Figure 10 suggests that the nonzero slope of this fit is a result of the two Hyades members with the lowest 
FUV fiux, i.e. the two points in the lower right corner of Figure 1101 Fitting the evolution while ignoring 
these two points gives 

logi?puy = (-0.33 ±0.16) log (age/yr) + (-1.95 ± 1.32) 
Thus, our fit is not significantly affected by the presence or absence of these two points. 



4. Testing for Systematic Errors 

The relations described by Equations and [MTUl and the R^jy indices require a host of assumptions. 
In this section, we verify that our results do not depend on these assumptions. In Section l4.1l we perform 
a simulation to test whether our procedure for fitting upper limits to the UV fluxes bias the resulting 
fits. In Section we explore whether changing the degrees of the polynomials in Equations [21151 and CTTUl 
significantly changes the colors, ^jjj^ indices, or ages predicted by the fits. Finally, in Section IT^ we explore 
whether our R'tj-^j indices depend greatly on the color of the source. 



4.1. Testing the Effects of FUV Incompleteness 



In Sections 12.31 and 13. 2[ we used an expectation-maximization (EM) algorithm presented by [Aitkin 



(|1981l ) to fit the UV flux upper limits in our data. The algorithm assumes that the scatter in the data 
is Gaussian, uniform across the sample, and intrinsic to the system under study rather than the result of 
measurement error. All three assumptions are violated in our data. Here we show, based on Monte Carlo 
simulations, that these assumptions do not invalidate our fits. 

We focused on reproducing Equation [2] and the left panel of Figure [21 as this is the result that seems the 
most suspicious: the UV-i?jjj^ relation levels out at the same B-V color at which we start missing stars in 
the FUV, suggesting the apparent fiattcning is an artifact of our detection limits. We therefore constructed 
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a model where the intrinsic UV-i?jjj^ relation remains linear in _B — V at all colors, based on a fit to stars 
with B — V < 0.7 where our observations are complete: 

{FUV - V)true = H-IQ - 2.17((logi?ijK)true + 4-5) - lA6{{\ogR'^^}true + 4-5)' 

+ 10.31((i? - y We - 0.8) - 7.26((S - Vk^ue - 0.8)((log i?^j^)true + 4.5) 
- 8.83((i3 - y)true " 0.8)((log i?^K)true + 4.5)^ (12) 

The linear B — V dependence of this model is in contrast to Equation [21 which has a quadratic dependence 
on B -V, leveling out over 0.7 < B - V < 0.9. 

We created a set of 1200 simulated stars by independently drawing (Vj^j-^g, {B — V^)true' (^^S-^HK-'true) 
triplets from the Hipparcos sample. Use of the Hipparcos data in place of analytical models ensured our 
model had similar distributions to our data. We then applied Equation [12] to get intrinsic FUV magnitudes 
FU\\Yue- 10% of the FUV magnitudes were brightened by an amount drawn from an exponential distribution 
with a mean of 1.5 mag. Introducing these bright outliers reproduced the "cloud" of UV-bright points seen 
in Figure[5] Aside from the population of outliers, we assumed there was no scatter in the intrinsic UV-i?jjj^ 
relation. 

To simulate the observing process, we randomly drew an exposure time and FUV sky background count 
rate for each star from the real GALEX observations. The exposure time and background were drawn 
independently from each other, and independently from the {V,B — VjlogR'^^) triplets, so that a typical 
simulated star was based on three different stars from the original data. Again, drawing model parameters 
from the data guaranteed the distributions we were using had the same properties - such as a tail of very long 
exposures - as the real data. We generated observed FUV fluxes for each star from a Poisson distribution 
with a mean equal to the intrinsic flux plus the background, multiplied by the exposure time. We calculated 
detection limits from Equation [1] and reported non-detections where the observed magnitude FUV^^^ was 
fainter than these limits. We generated observed V magnitudes, B — V colors, and logR'^^ measurements 
simply by drawing Gaussian random numbers with means V^^uei ~ ^)true' ^^'^ (^^g i?jj-j^)^]-yg) and 
standard deviations 0.01 mag, 0.01 mag, and 0.08 dex respectively. We then fit polynomials to the data 
following the procedure of Section 12. 3[ including the use of the leave-onc-out cross-validation score as a 
goodness-of-fit statistic to determine the degree of the polynomial. 

The best fit by cross-validation score is linear in B — V and quadratic in log ^jjj^i like Equation [12] The 
best fit that allows curvature in i? — y is quadratic in _B — and linear in log -Rjjj^ j shown in Figure 1111 
Figure [TT] looks very different from Figure [5] the equivalent plot in our real data. We believe the lack of 
curvature in the simulated fit arises because our fits to the simulated data are very poorly constrained at red 
B — V colors. Our simulated FUV observations detect only 37 of 394 stars (9%) with B — V > 0.8, and 6 of 
177 stars (3%) with B -V > 0.9, compared to our real detection rate of 81 of 396 (20%) with B -V >0.8 
and 29 of 184 (16%) with B — V > 0.9. This is because Equation [T2l places most cool stars 2 magnitudes 
below our sensitivity limits, so only the "outliers" - the 10% of stars that were artificially given an FUV 
excess - are detected. To get the simulations to reproduce the observed detection rate at red B — V, we need 
to artificially brighten ~ 60% of the FUV magnitudes rather than 10%. Not only does brightening so many 
stars increase the number of UV-bright outliers at blue B ^ V colors well beyond the observed number, it 
means that Equation 1121 no longer describes the simulated population as intended. 

We conclude that the intrinsic FUV — vs. B — V relation does flatten at B — V > 0.8, as otherwise 
we would expect fewer detections and a qualitatively different fit than we observe. However, since Figure [TT] 
shows the fit to the simulated data is shallower than Equation [121 it is still possible that Equation [2] 
exaggerates the flatness of the intrinsic FUV — V vs. B — V relation. 



- 14 - 



4.2. Testing the Dependence on Functional Form 
In Sections 12.31 and 13.21 we fit our data witli bivariate polynomial s, choosing the de gree of the polyno- 



mial that minimized the leave-one-out cross validation score CV (e.g. I WassermanI 120061 ) . which acted as a 
goodness-of-fit statistic. However, the solution with the lowest CV score was often one of several with very 
similar goodness of fit. Since CV, because of its generality, does not follow a specific distribution such as 
the chi-squared distribution, we cannot judge whether a difference in CV is statistically significant. This 
potentially makes the choice of the best form for the fit sensitive to fluctuations in CV. 

Fortunately, the functional form of the fit - in our case, the degree of the polynomial - carries no physical 
significance. The decision to choose one form over another is only important insofar as it changes the value 
of UV — V, R'hk^ NUV — J, or age predicted for a particular star. In this section, we show that replacing 
the polynomial fit with the lowest CV score with, for example, the fit that gives the second-lowest score does 
not significantly change the numerical values predicted by the fit. It follows that fiuctuations in CV do not 
matter: all the fits with similar CV scores make the same predictions. 

For each of Equations [2][S] and IWlUl we compared the fits we presented in Sections 12.31 and 13.21 which 
had the lowest CV score, with the alternative polynomials that had the second, third, and fourth lowest 
scores. For example, our Equation [5] presents FUV —V as a. quadratic function oi B ~V and a quadratic 
function of logi?jjj^. The model that had the second- lowest CV score gave FUV — as a cubic function 
of -B — and a quadratic function of logi?jjj^, but with best- fit parameters that caused the two models to 
predict the same value of FUV — V, to within 0.03 mag, over most of the B — V vs. log R'^iK P^^^'^^'^ter 
space. It follows that the choice of polynomial form for Equation [51 quadratic or cubic in B — V, matters 
only at the ^ 0.03 mag level. 

In Table [21 we present the RMS difference between each model adopted in Equations [2l[5l and iQlfTOl and 
each of its three alternatives, with the RMS taken over the region in parameter space where each model 
is valid. For comparison, we also list the residuals we found around the equations. The RMS differences 
between alternate forms for the fits are typically a factor of 3-4 less than the scatter around the fits, indicating 
that choosing one of the alternate forms instead would not have changed our results. 



4.3. Testing for Selection Effects 

The trend between UV color and -Rjjj^ seen in Figures [2E] is small, comparable to the scatter in the 
data. Here we show that this trend is not the result of flux biases, such as active stars being more distant 
or less luminous and therefore having UV fluxes more prone to measurement error. 

We binned the data in logi?jjj^ and B — V and calculated the average apparent V magnitude and 
average parallax for each bin. We present the means and the errors on the means in Table [31 In each B — V 
range, we used a two-sided t test to test whether the means in different log ^jjj^ bins were significantly 
different. No pair of mean parallaxes differs at more than 82% confidence. For stars with B ~ V > 0.7, 
only one pair of mean magnitudes differs at more than 95% confidence, which is not a significant result in 
9 tests. However, for stars with B — V < 0.7 the least active bin (log-Rjjj^ < —5.0) differs from either of 
the other two bins at 99.9% confidence or greater. Inactive stars with 0.5 < B — V < 0.6 are on average 0.5 
magnitudes brighter than their active counterparts, while those with 0.6 < B — V < 0.7 are 0.3 magnitudes 
brighter. This correlation represents the fact that evolved main sequence stars tend to be both less active 
and more luminous than their ZAMS counterparts. 
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Since the only correlation between apparent magnitude and -Rjjpj is at the blue end of our sample, 
whereas our inferred relationships between UV flux and ^jjj^ are strongest for the reddest stars in our 
sample, we conclude our results are not caused by a bias in stellar magnitude with -Rjjj^- 



4.4. Testing for the EfTects of Source Color 

In Section 12.3.31 we inferred excess UV fluxes from our GALEX photometry. The inferred fluxes are 
proportional to the instrument count rate CPS — ^ Af^^[v)fu{v)/{hv)dv, where is the effective area 
for the FUV or NUV detector. However, we interpret the fluxes as the amount of emission in a well-defined 
band, /uy — fu{v)di', which does not correlate perfectly with count rate. Here we show that the color 
of the source does not introduce large inaccuracies in our inferred /uV- 

We considered three t emplates for f, /. a 10 "^ K blackbody, a 10^ K blackbody, and a list of FUV lines 



in the a Gen spectrum of iPagano et al.l (|2004l ). The first two models bracket the temperature range at 
which most FUV and NUV emission lines and continuum form. The third allows us to consider a third 
extreme, a spectrum dominated entirely by lines with no continuum contribution. Stellar UV spectra are a 
superposition of 10^~^ K emission from both lines and continuum, so analyzing the color corrections for the 
three models places an upper limit on the variation we expect to see in real data. 

Evaluating the integrals numerically, we normalized each template spectrum to give the same GPS in 
the FUV or NUV band, as appropriate. We then evaluated /uy directly from the scaled spectrum. We found 
that the 10^ K blackbody that produces the same FUV GALEX response as the lO** K blackbody has 92% 
the latter's flux in the 1350-1780A range. The a Gen line list has 97% the flux of the 10"^ K blackbody. 
The 10^ K blackbody that produces the same NUV response as the 10"' K blackbody has 107% the flux 
in 1780-2830A. We could not evaluate the NUV response to the a Gen line list because it only extends to 
1690A. 

In summary, the true flux associated with a GALEX magnitude varies by up to 8% over the range of 
templates we have adopted. Since the three template spectra bracket the variety of source colors we expect 
to encounter, we likewise expect that color systematics will not cause the true flux to deviate by more than 
8% from our estimate. 



5. Discussion 

5.1. UV Flux Evolution v^rith Age 

In our data, both GALEX bands probe multiple physical environments. The NUV flux is dominated by 
the photosphere, as seen in the close agreement of synthetic and observed NUV colors in Figure [5l However, 
the stellar chromosphere, i.e. optically thin gas cooler than ^ 1 — 2 x 10'* K, also contributes to the NUV flux 
through lines from a variety of low-ionization species. The FUV flux is, to order of magnitude, evenly split 
(see Appendix |B]) between continuum emission from the chromosphere and lines from both the chromosphere 
and the transition region (a narrow zone of lO'* K to ^ 5 x 10^ K gas). Our results on stellar activity in the 
UV need to be interpreted with these mixtures in mind. 

The most active stars in our Hipparcos sample typic ally show -Rp-gy "^. -^ .Q""^'/ those in the moving 
group and cluster sample can reach up to -Rpuy ~ lO""^'^. iGorti &: Hollenbacbl (|2009l ) found, by integrating 
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over spectra from Valenti at al. ( 2000[) . that weak-lined T Tauri stars show a median R'g^^ 2070A ^ 
Their figure quahtatively agrees with our results, if we assume that 50-80% of the excess flux in 913-2070A 
is emitted by the Lyman a line, and if we also note that our FUV band (1350-1780A) is about one third as 
wide as theirs. 

Our relations between the FUV excess flux and the Ca II H and K line flux are likewise consistent with 
previous work on co rrelations between the Ca II and UV lines. The most extensive such comparison was by 
Rutten et al.l (|l99l[ ). who found: 



logi^Si II 
log Fen 
log^si IV 
logFc IV 



1.36 log Fhi^ - 3.8 
1.65 log Fhi^ -6.2 
2.141ogFiiK -9.3 
1.851ogFHK-7.2 



where F denotes the flux in a particular line at the stellar surface. iRutten et al.l did not subtract a contribu- 
tion from the photosphere when calculating their Ca II H and K surface fluxes, so to compare to their results 
we also avoid subtracting a photosphere contribution. Following the procedure we used to fit Equation [8l 
but using the index Rn}^ to include the photosphere fiux, we find 



log R'p^y = (1.87 ± 0.07) log%K + (3-30 ± 0.33) 



(13) 



This relation has a similar slope to the slopes iRutten et al.l (|l99lh found for four prominent transition region 
lines, even though only part of the FUV flux originates in the transition region. 

Our work in the 1350- 1780 A range fills a gap in previous studies of the evolution of stellar activity with 



age. iRibas et al.l (|2005lO presented power-law fits to the decay of X-ray and ultraviolet activity with age 



among solar analogs. They noted that between the ages of 100 Myr and 6.7 Gyr relatively long-wavelength 
emission falls off much more slowly with age than shorter-wavelength emission: soft X-ray (0.6-12 keV) 
emission falls off as t'^-^'^, extreme UV (100-360A) emission falls off as r-^ ^", and far-UV (920-1180A) 
emission falls off as r~° *^ Our Equation [TT] implies that between 30 Myr and 1 Gyr 1350-1780A emission 
falls off as r^^-'^'^^'^-^^ , continuing the trend although with large scatter. We caution, however, that because 
we could not estimate R^UY stars of cooler spectral type than K5 this result is based primarily on the 
AB Dor moving group and the Hyades cluster, as these groups contributed the bulk of the G- and K-type 
stars in our sample. We also note that Equation [8] shows that the excess FUV and excess Ca II fluxes are 
proportional to each other, suggesting that the trend of long-wavelength emission decaying more gradually 
does not extend to longer wavelengths than the FUV. 



5.2. Predictions Among Rj^j^i Stellar Age, and UV Color Excess 

5.2.1. The Residuals and their Significance 

The large residuals around our relations between the Ca II activity indicator ^in^, stellar age, and UV 
color excess. Equations [2]|5l and [MTOl reflect a genuine limitation of the data. Stars with very different Ca II 
fluxes can show similar UV fluxes, and vice versa, as suggested by Figures [3]|4] and [SEl To show that the 
scatter is not an artifact of our fltting methods, or even of the decision to fit curves in the first place, we 



- 17 - 



select from the Hipparcos sample pairs of stars whose colors are identical within la confidence: 



{{FUV - V)2 - {FUV - V)iY 



FUV,1 



FUV,2 ^"V,2 



{{B-Vh-iB^V),r 



< 1 



{{NUV - V)2 - {NUV - V)if 



NUV,1 ^YA 



'NUV,2 "^"V,2 



aB-Vh-iB-V),Y 
'^B-Y,l + '^B-V,2 



and 



< 1 



(14) 



where a denotes the formal uncertainty on the color or magnitude, and the subscripts 1 and 2 refer to the 
two stars in each pair. Our Hipparcos sample of 1204 stars allows for 724,206 possible pairings. Of these, 
1875 meet the criterion in Equation [TH We will refer to these 1875 pairs of stars as photometric analog 
pairs. Some stars are members of multiple analog pairs, i.e. their colors resemble those of more than one 
other star. Analog pairs typically differ in B — V color by a median of 0.017 mag, in FUV — V color by 
0.14 mag, and in NUV — V color by 0.016 mag. 

Since the two stars in each analog pair have statistically indistinguishable photometry, one might hope 
that they also have similar Rjj^ measurements. We show the actual distribution of log -Rjjj^ 2 ~ ^'^S ^HK 1 
for all 1875 pairs in Figure [T^ A quarter of the pairs indeed have identical logR'^^ values to within our 
observed scatter of 0.08 dex, but another quarter of the pairs have highly discrepant logi?jjj^ values differing 
by ~ 0.3-0.5 dex. The median difference in logi?jjj^ is 0.17 dex, close to our RMS residuals of 0.18 dex 
around Equation \5\ We have not performed a fit of any kind to arrive at this number; we have simply 
picked out photometrically identical stars from the data and asked whether their ^jjj^ measurements are 
also consistent with each other. 



The predictive power of all six of our fits - Equations [2||5] and IMTOl ~ is limited by these large residuals. 
Equations [21 HI and [9l which predict the UV fiuxes of stars based on their Rjjj^ values and optical or 
near-infrared colors, have residuals of 0.33 mag, 0.21 mag, and 0.46 mag, respectively. If we assume that 
FUV — V for all stars with similar ^jjj^ follows a Gaussian distribution with a standard deviation equal 
to our 0.33 mag residuals, we would report a 95% (2cr) confidence interval for a star's FUV magnitude that 
allows nearly a factor of 2 variation in fiux. These predictions, uncertain as they are, are still much more 
precise than those we would get were we to ignore stellar activity. If we repeat the fitting procedure for 
Equations [21 m and[9l but require that the models have no dependence on ^^jjj^ or age, we find much larger 
residuals of 0.68 mag, 0.25 mag, and 0.64 mag, respectively. 

While Equations [21 [4l and [9| allow factor of 2 predictions in UV flux, the residuals around Equations [3l 
[SI and [TU] (0.15 dex, 0.18 dex, and 0.39 dex, respectively) prevent us from using UV fluxes to predict 
^HK better precision than half its observed range, or stellar age to better than an order of magnitude. 
Consider as an example a star with B — V = 0.65 and FUV — V = 11.25, near the UV-bright edge of 
the stellar locus in Figure [H Equation [3] states that such stars have log Rjjj^ = —4.48 on average. If 
we assume that log ^jjj^ for all stars with similar colors follows a Gaussian distribution with a standard 
deviation equal to our RMS residuals of 0.15 dex, we report a 95% confldence interval of —4.78 < log ^jjj^ < 
—4.19. This is an enormous ran ge in activity, corresponding to an age anywhere between 50 Myr and 3 Gyr 
( Mamaiek &: Hillenbrand 2008[l . Even though the star has a very high UV flux compared to other stars of 
the same B —V , we can infer little more than that the star is more active than the Sun. We found a similar 
result in ^ Findeisen fc Hillenbrand ( 2010 ). where we selected active stars by their UV flux but found that 
only a third showed obvious Ca II or Ha emission in follow-up optical spectra. The remainder, like the 3 Gyr 
stars in the preceding example, were not unusually active. 
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5.2.2. Systematics in Our Fits 

Equations [5][5] can give counterintuitive results if used carelessly, because they are best-fit models in 
the nonparametric sense. Instead of specifying a hypothetical one-to-one relationship between UV flux and 
log ^jjj^ , they give the average UV — V in a large sample of stars with the same measured values oi B — V 
and ^jjj^j or the average logi?jjj^ in a large sample of stars with the same measured values of i? — and 
UV — V. The average is, implicitly, over a population with the same ^jjj^ distribution as our sample. The 
log-age predictions of Equations CTTUl must be treated as averages in the same sense. This averaging process 
is unusually significant for our fits because they have large residuals, and most important for Equations [3J 
m and llOi where the residuals are a significant fraction of the width of the parameter space. 

One consequence of the fits acting as averages is visible in Figure El Equation [10] systematically under- 
predicts the age of the oldest stars in the moving group and cluster sample, those in the Hyades, and 
systematically over-predicts the age of the youngest stars, those in the TW Hydrae association. The reason 
can be seen with a thought experiment. A particular (high) observed UV flux can be produced by a very 
young star, such as a TW Hydrae member, with UV emission typical for its age, or it can be produced by 
a moderately young star, such as a Tucana/Horologium or AB Doradus member, that is (or is measured to 
be) UV-bright for its age. It cannot be produced by an even younger star that appears UV-faint for its age, 
because TW Hydrae is the youngest association in our sample. The mean age predicted by Equation 1101 for 
a star with the UV flux typical of a TW Hya member will be older than the age of TW Hya, because there 
are older stars with the same flux but no younger stars. More generally, the ages predicted by Equation 1101 
are biased away from the edges of the 10^'^ — lO^'* Myr age range probed by our moving group and cluster 
sample. 

The edge effect we identifled for ages is much weaker for the -Rjjj^ values predicted by Equations [3] and 
m because the Hipparcos sample is volume limited so the activity distribution in the sample approximates 
that in the Galactic stellar population. However, the most extremely active (log-Rjjj^ > —4.0) or inactive 
(logi?jjj^ < —5.5) stars are not represented, so we expect Equations [3] and [S] will slightly underpredict the 
^HK '^^ most active stars and overpredict that of the least active stars. 

Other systematic effects may also skew the predictions at very hi gh or very low activit y levels. We tested 



the predictive power of Equations [3] and [5] on the activity sample of iDuncan et al.l (|l99ll ). They gave Ca II 
fluxes, in the form of the Mount Wilson 5-index, for 12 stars that were detected by GALEX but were too 
faint to include in the Hipparcos catalog. Because they are not Hipparcos stars, these 12 stars form a sample 
independent of our calibration data. 8 of these stars meet the requirements of either Equation [3] or [5l We 
found 95% confldence intervals for logi?jjj^ based on the GALEX data, approximating the scatter around 
our flts as Gaussian in logi?jjj^. In Table El we compare our confidence interval s with log Rjjy va lues we 



computed from the published S'-index measurements following the prescription of lNoves et al.r(|1984l ). Only 
5 of the 8 stars' log_Rjjj^ measurements lie within our 95% confidence interval, while we expect at least 7 
out of 8 to do so. Since the scatter around the fits is close to Gaussian, as assumed, it follows that there 
are as yet unidentified systematics in the fits themselves. We note that the three discrepant stars include 
the two most active stars in the sample, with logi?jjj^ > —4.0, suggesting that very active stars may be the 
problem area. 
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5.3. Unaccounted-For Variables 



In this study, we have ignored nietalhcity. We expect metaUicity to affect stellar UV luminosities in three 
ways. First, we expect the photospheres of metal-rich stars to have highe r UV opacit ies, particularly from line 
blanketing, and redder UV-optical colors compared to metal-poor stars (jCravlboOSh . Second, FUV emission 



comes from metal lines and continuum, so chromospheric and transition region gas around more metal-rich 
stars should produce more FUV flux per unit emission measure. Finally, metal line emission is the primary 
cooling mechanism for transition region gas, so we expect the structure of the outer stellar atmosphere - 
in particular, the thickness and density of the transition region - to depend on the stellar metallicity, with 
potentially complex effects on the observed activity. The number of processes to untangle, and our inability 
to thoroughly sample a large metallicity range, preclude our study of the metallicity-dependence of stellar 
UV flux. 

We have also ignored rotation, even though activity-rotation relations are, on theoretical grounds, more 
fundamental than activity-age relations. Stellar activity results from chromospheric and coronal heating, 
which among Sun-like stars is mediated by the s tellar magnetic field . Every major model for the solar dynamo 
incorporates differential rotation in some form ( CharbonneaTjIiolO ) . and so the strength of the magnetic field 



might be expected to scale with the overall rotation rate of the star. Activity should therefore be strongly 
correlated with rotation, whereas the relationship between activity and age represe nts the combination of the 



activit y-rotation correlation with the systematic slowing of stellar rotation with age. iMamaiek fc Hillenbrand 



(|2008[ ) confirm that activity-rotation-age relations show a tighter correlation than activity-age relations that 
do not incorporate measured rotation periods. However, rotation periods ar e not available for as many stars 
as ^jjj^, and are not even detectable in all stars surveyed for rotation (e.g. Donahue et al. 19961) . 



While we expect stellar UV fiux to be affected by both metallicity and rotation, we do not need to 
invoke either factor to explain our data. The scatter we observe in UV flux is consistent with that expected 
from propagated errors and previously measured variability, so without more precise photometry and coeval 
activity measurements the data do not support adding more variables to the analysis. We have instead 
restricted the problem to its three essential dimensions: UV flux, photosphere effective temperature, and 
one age or activity indicator to which we may compare the UV flux. 



6. Summary 

We set out to quantify the dependence of stellar ultraviolet broad-band flux on age among Sun-like 
stars, and to test the correlation of UV flux with more traditional activity indicators, in particular the well 
studied chromospheric indicator -Rjjj^- To this end, we matched both an unbiased, volume-limited sample 
of nearby stars with ^jjj^ measurements but unknown ages, and an age-calibrated but incomplete sample 
of young stars, to archived GALEX data. To avoid introducing a model dependence into our results, we 
used nonparametric fits to describe the correlation between UV-optical or UV-infrared colors and ^jjj^ or 
age measurements. We also explored the construction of model-dependent activity indicators, ^p^V ^^'^ 
'^NUV' '^'^'^^^'^ analogously to ^jjj^- 

Our major results may be summarized as follows: 

• We detect a clear correlation between UV color and both ^jjj^ and age, in the sense that younger 
stars with higher -Rjjj^ have higher UV luminosities. We provide relations describing the average -Rjjj^ 
of stars at fixed UV color, and the average UV color of stars at fixed Rjj^, where both averages are 
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understood to be taken over a volume-limited sample. We likewise construct relations describing the 
average age of stars at fixed UV color, and the average UV color of stars at fixed age, although given 
the heterogeneity of the sample on which these latter relations are based the meaning of the averages 
is not as clear. 

Our mean relations are not good predictors of individual ^jjj^ values (±0.15 — 0.18 dex) or ages 
(±0.4 dex, with additional systematics). Selecting young, active stars using UV excess techniques 
requires a thorough s pectroscopic follow-up campaig n that is likely to reject many UV-selected targets 
as false positives (cf. Findeisen fc Hillenbrandl 2010l ). Variability and measurement error appear to be 



equally responsible for the scatter around our mean relations. 

We find that attempts to define ^^puy or -Rj^yy indices are dominated by systematic uncertainties 
in the photosphere contribution, particularly in the NUV. Until these uncertainties can be constrained 
observationally, use of such indices is likely premature. We recommend instead a purely empirical 
activity index such as UV-optical color, despite the strong temperature-dependence such indices have. 

We find that the FUV excess fiux between 1350A and 1780A is proportional to the excess flux in the 
Ca II H and K lines. We tentative l y find that the FUV excess flux decays with time as t--0-42±o.18^ 



extending the results of lRibas et al.l (j2005[ ) to longer wavelengths. This decay rate should be conflrmed 



with data that sample a wider variety of ages. 

Our results on the age dependence of the FUV flux are currently limited by the lack of young stars of 
known age with GALEX data. Improving these results will likely require expanded moving group membership 
lists - most open clusters near enough to probe down to K-type photospheres at AIS survey depth have at 
least one star bright enough to pose a danger to GALEX, forcing GALEX to either observe only the outskirts 
of the cluster or avoid it entirely. Acquiring more precise near-infrared photometry than the often saturated 
2MASS fluxes we used would also improve many of our age results, even if the stellar sample remained 
unchanged. 

^The GALEX Afl Sky Survey contains 111,755,312 sourcetH, of which 80% are stars ( Bianchi et al 



20051 ). Far more stars already have ultraviolet photometry than will have spectroscopic measurements in the 
foreseeable future. We have shown that ultraviolet photometry does not show a tight correlation to optical 
activity indicators such as ^jjj^i and cannot be used to cleanly identify high-i?jjpr stars. However, it can 
still be used to filter a sample, removing the most inactive stars to allow more efficient follow-up of the rest. 
No other activity indicator is as readily available as FUV and NUV photometry. 



We would like to thank the referee for many insightful comments and suggestions, Ted Wyder and Patrick 
Morrissey for addressing our questions about the GALEX data, and Chad Schafcr and John Carpenter for 
their advice on the statistical analysis. This research was supported in part by NASA grant NNX08AH95G 
to L.A.H. 
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A. Optimal Photosphere Colors for UV Studies 



As we noted in Findeisen &: Hillenbrand ( 2010f) . the precision to which UV photometry can be used as 



an activity measure is limited by the precision to which we know the photosphere UV flux. The uncertainty 
in the photosphere flux, in turn, is dominated by the precision to which the star can be classified because 
of the steep dependence of UV flux on spectral type. We consider here the use of a single color for stellar 
classification, implicit in the use of color-color diagrams such as Figure [2l This analysis does not apply to the 
multi-band fits we used to derive R'jjy ™ Section 12.3.31 However, it would apply to an attempt to use the 
B — V color to estimate both m-^Qi and the photospheric UV fiux, and derive Rjjy from these single-color 
values. 



The uncertainty in the UV flux is roughly 



(d{UV-X) 



"UV - — Qc — 



/ AjUV^X] 
\ AC 



crc 



= AiUV-X)^ (Al) 

where the color UV — X represents the UV fiux relative to flux in a purely photospheric band, C is the 
color used to classify the star, and A denotes the range of a quantity observed along the main sequence. 
Minimizing the first factor in Eg nation I Al I involves choosing a diagnostic color UV — X that varies relatively 
little with spectral type, while minimizing the second involves choosing a photosphere color C that resolves 
the main sequence into as many elements as possible. 

A{UV — X) is smallest if X is as blue a band as possible. In particular, an analysis based on UV — B 
or UV ~ V color will always be more precise than one based on UV — J or UV — K. We adopted UV — V 
rather than UV — B because we could determine stars' V magnitudes much more precisely, and because the 
penalty Equation lAl I predicts for using V, '^[piryZ^] a[nuv~^] ^-^^ small. 

Minimizing the classification term ac/ AC is more complex, because it depends on the quality of the 
available photometry and not just on the position of the stellar locus in color space. We found ac and AC 
for several colors in our Hipparcos sample, adopting the median formal error for a and adopting the range 
between the 10th and 90th percentiles for A to avoid biases from outliers. We found as-v / ^{B — V) =0.04, 
uj-kI A{J — K) = 0.35, and av-K/A{V — K) — 0.08. Infrared colors are a poor classifier for our sample 
because many of our target stars are saturated in 2MASS, giving large (~ 0.2 mag) errors in their near 
infrared magnitudes. As a result, we identified stars by their B — V color for the Hipparcos sample, even 
though using J — K would allow us to compa re directly to our results for th e moving group and cluster 
sample, or to the measurements we obtained in iFindeisen fc Hillenbrandl (|2O10l ). 



B. Expected Contributions to the FUV Flux 



Pagano et al.l (j2004l ) present a detailed STIS spectrum of a Cen A, including a complete line list (their 



Table 4) and continuum estimates (dashed line in their Figure 2). We use this spectrum as a template to 
estimate the contribution to the FUV flux from various sources in Sun-like stars. 

Adding together the fluxes of all the lines between 1350A and 1690A, we find a Cen A has a total FUV 
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line flus0 of 2.3 x 10"" erg s"^ cm-^ 



5 X 10 ere 



erg s " cm " . The continuum seen in iPagano et alj (j2004 ) Figure 2 varies between 
rg s" " cm'^ at 1370A and 3 x 10"^"^ erg s~^ cm~^ at 1690A. Assuming an average 
continuum of 1.5 x 10~^^ erg s~^ cm~^ A~^, we infer a total continuum flux of 5.1 x 10~^^ erg s~^ cm~^ 
over the range 1350-1690A. It follows that, for a typical Sun-like star, lines contribute 30% of the GALEX 
FUV flux, with the rest coming from a weak continuum. It is possible that even at the high resolution of 
the STIS observations (_R ~ 114, 000) some of the "continuum" consists of unresolved lines, but this has not 
been investigated (Pagano, priv. comm.). 

Although the Lyman a line in thei r spect rum is heavily contaminated with both interstellar absorption 



and geocoronal emission, iPagano et al.l (j2004l ) estimated a flux of 1.04 x 10 erg s ^ cm ^ by fitting the 



wings of the line profile. In comparison, the total line flux over the range covered by the observations, 
1170-1690A, is 1.39 x 10~^° erg s~^ cm^^, and our estimate for the continuum in that range is 7.8 x 
10~^^ erg s~^ cm~^. Therefore, Lyman a. contributes 75% of the line flux and 48% of the total flux in the 
range 1170-1690A. 
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Fig. 1. — UV vs. B — V diagrams of our volume-limited Hipparcos sample, showing the effect of satura- 
tion. The dotted line shows the magnitude at which nonlinearity effects introduce a 10% or larger error 
(jMorrissev et al.l 120071) . The ellipses in the upper right of the left panel and the lower left of the middle 
panel show the median errors in flux and color. The two right panels show the stars that did (center) and 
did not (right) pass the requirement that V + 6A6{B — V) > 12.8. This cut effectively removed stars that 
saturated in the NUV, without introducing a bias toward stars that were underluminous in the NUV. Our 
FUV data does not saturate, so no cut was necessary. 




Fig. 2. — UV — y vs. B — V diagrams of our volume-limited Hipparcos sample. Red points are the most 
active stars, with logR'^j^ > —4.5. Green points have —5.0 < logR'^j^ < —4.5, while blue points have 
logi?jjj^ < —5.0. The ellipses in the upper left corner of either panel show the median errors in color. To 
illustrate the dependence on activity, we plot Equations ^ and U at the median log R'ttt^ of each bin. 
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Fig. 3. — UV — V vs. ^jjj^ diagrams of our volume-limited Hipparcos sample. In the left panel, red points 
have B — V > 0.8, green points have 0.65 < B ~ V < 0.8, and blue points have B — V < 0.65. In the right 
panel, red points have B — V > 0.9, green points have 0.8 < B — V < 0.9, and blue points have B — V < 0.8. 
The ellipses at the top of either panel show the median errors in color. The utility of the UV as an activity 
indicator can be seen directly as a downward trend in UV — V color with activity. We plot Equations [5] and 
|4] at the median _B — of each bin. 
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Fig. 4. — ^jjj^ "^s- UV — V diagrams of our volume-limited Hipparcos sample. The left panel only shows 
stars with V + 12.03{B — V) < 15.5, where our FUV observations were complete. In the left panel, red points 
have B ~ V > 0.7, green points have 0.6 < B — V < 0.7, and blue points have B — V < 0.6. In the right 
panel, red points have B ~V > 0.9, green points have 0.8 < B — V < 0.9, and blue points have B — V < 0.8. 
The ellipses at the top right of either panel show the median errors along each axis. The main difference 
between these figures and Figure [3] is that we plot Equations [3] and [SJ which give the average -Rjjj^ at fixed 
UV — V rather than the average UV — V at fixed R'ttt^ ■ 
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Fig. 5. — Color-color diagrams comparing our synthetic UV fluxes to our observations with GALEX. Colors 
have the same meaning as in Figure [2l The dotted curve shows the stellar locus predicted by matching the 
synthetic FUV and NUV rn agnitudes calculated for each model temperature with the photometry listed by 



Kraus fc Hillenbrand! (j2007l) for the same temperature. The results differ markedly from the observations, 



particularly in the NUV, where the synthetic magnitudes imply that a large fraction of stars emit less UV 
flux than their photospheres can account for. The solid curve shows the stellar locus after the photosphere 
FUV and NUV magnitudes have been corrected for temperature systematics as described in the text; the 
results appear much more consistent with the NUV observations. 
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Fig. 6. — -^uv ^HK P^*^^® '^^ °™ volume-limited Hipparcos sample. The left panel, a log plot, shows 
the correlation between R'^uy ^^i^ ^HK' "^^"^ solid line represents the fit given by Equation |8l The right 
panel, a linear plot, shows that ^p^jjy measurements are very noisy with many negative values. The ellipse 
in both panels shows the median errors. The distribution of ^^TyrTj-ir is centered on zero. 
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Fig. 7. — UV — J vs. J — K diagrams of our cluster and moving group sample. Each symbol corresponds 
to a different association, as shown in the legend, while the hue maps linearly to log age. The ellipses at the 
upper left of either figure show the median errors in color. The NUV panel shows a clear trend with age, 
but any relation between FUV and age is obscured by non-detections. 
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Fig. 8. — UV — J vs. age diagrams of our cluster and moving group sample. Symbols are the same as in 
Figure[71 In the left panel, red points have J — K > 0.7, green points have 0.4 < J — K < 0.7, and blue points 
have J — K < 0.4. In the right panel, red points have J — K > 0.8, green points have 0.6 < J — K < 0.8, 
and blue points have J — K < 0.6. To illustrate the slight trend with age, we plot Equation[9]at the median 
J — K oi each bin. 
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Fig. 9. — Age vs. UV — J diagrams of our cluster and moving group sample. Symbols are the same as in 
Figure[7l The left panel only shows stars with J + 17.29(J— K) < 13.87, where our FUV observations were 
complete. In the left panel, red points have J—K > 0.3, green points have 0.2 < J~K < 0.3, and blue points 
have J — K < 0.2. In the right panel, red points have J — K > 0.8, green points have 0.6 < J — K < 0.8, 
and blue points have J — K < 0.6. The main difference between these figures and Figure [8] is that we plot 
Equation [TOl which gives the average age at fixed NUV — J rather than the average NUV — J at fixed age. 
While UV-bright stars do tend to be older than UV-faint stars on average, there is so much scatter at fixed 
age that the age of any individual star cannot be predicted well. 
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Fig. 10. — vs. age plots of our moving group and cluster sample. Symbols are the same as in Figure [T] 

The left panel, a log plot, shows the correlation between ^pyy and age. The solid line represents the 
fit given by Equation [TT] The right panel, a linear plot, shows that -Rj^yy measurements appear to be 
dominated by noise. The filled circle in either panel is not a real point, but illustrates the median error in 
RjjY' Both panels include only stars whose inferred spectral type is K5 or hotter, as we lack the models to 
estimate photospheric UV fluxes for cooler stars. Thus, there are fewer points than in Figure |51 in particular 
none younger than ~ 30 Myr. 
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Fig. 11. — A typical run of our Monte Carlo simulations in Section HTTI The panel on the left shows the 
intrinsic color-color relation, without observational errors or selection effects, together with Equation [T2l 
which wc used to generate the data. The panel on the right shows the same sample, but with observational 
errors and flux limits applied, and a fit to the simulated data. On both panels red points are the most 
active stars, with log-Rjjj^ > —4.5, green points have —5.0 < logR'^j^ < —4.5, while blue points have 
logi?jjj^ < —5.0. Equation [T^ and our fit to the simulated data are plotted at the median logi?jjj^ of each 
bin. 
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Fig. 12. — The fraction of pairs of photometric analogs differing by more than the amount along the abscissa 
in A log -Rjjj^- Although 26% have identical log ^jjj^ values within the typical scatter of 0.08 dex, 26% differ 
by more than 0.3 dex, and 16% differ by more than 0.4 dex. Because such large differences in log -Rjjj^ aic so 
common among photometrically indistinguishable stars, any attempt to predict log ^jjj^ from photometry 
alone will show large residuals, no matter what algorithm is used to make the predictions. 
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Table 1. Magnitudes Used to Fit Photosphere Fluxes to Observed Stars 



Spectral Type 


MpuY 




Mb 


My 


A/j 


AfH 


a/k 


^^bol 


B8 


0.80 


0.57 


-0.31 


-0.20 





01 





10 





11 


-1.00 


AO 


3.08 


2.26 


0.59 


0.61 





54 





58 





56 


0.30 


A2 


5.05 


3.30 


1.36 


1.31 


1 


12 


1 


15 


1 


12 


1.10 


A5 


7.42 


4.31 


2.06 


1.91 


1 


53 


1 


52 


1 


48 


1.75 


A7 


8.44 


4.83 


2.41 


2.21 


1 


75 


1 


71 


1 


66 


2.08 


FO 


10.50 


5.83 


3.02 


2.71 


2 


10 


2 


01 


1 


96 


2.61 


F2 


11.27 


6.29 


3.35 


3.01 


2 


32 


2 


20 


2 


14 


2.89 


F5 


13.41 


7.46 


4.17 


3.76 


2 


85 


2 


67 


2 


61 


3.61 


F8 


16.35 


8.77 


4.93 


4.41 


3 


31 


3 


08 


3 


01 


4.24 


GO 


17.65 


9.40 


5.24 


4.66 


3 


53 


3 


27 


3 


20 


4.47 


G2 


18.65 


9.88 


5.44 


4.81 


3 


64 


3 


38 


3 


30 


4.60 


G5 


19.74 


10.57 


5.80 


5.11 


3 


86 


3 


56 


3 


48 


4.89 


G8 


21.11 


11.56 


6.45 


5.71 


4 


31 


3 


95 


3 


86 


5.30 


KO 


22.66 


12.53 


6.83 


6.01 


4 


49 


4 


10 


4 


00 


5.69 


K2 


24.73 


14.02 


7.43 


6.51 


4 


80 


4 


35 


4 


24 


6.08 


K4 


28.23 


16.26 


8.16 


7.11 


5 


08 


4 


56 


4 


43 


6.55 


K5 


30.35 


17.72 


8.56 


7.41 


5 


20 


4 


64 


4 


51 


6.68 



Note. — FUV and NUV magnitudes are calculated on the GALEX system 
as described in the text. J, H, K, and bolometric magnitudes are taken from 



Kraus fc Hillenbrandl ( 2007 ) Table 5. B and V magnitudes are derived from the g 



an d r magnitudes in iKraus fc Hillenbrandl (|2007f ) Table 5 using the transformations 



of Jester et al 



(|2005l) . FUV and NUV magnitudes are not available for K7 or cooler 
spectral types because Kurucz models do not cover stars with redder B — V. 
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Table 2. RMS Change in Our Fits Under Alternate Functional Forms 



Adopted Fit 
(Model 1) 


Residuals 
(Model 1) 


(A/)rms 
Model 2 - Model 1 


(A/)rms 
Model 3 - Model 1 


(A/)rms 
Model 4 - Model 1 


Equation [2] 


0.33 mag 


0.033 mag 


0.018 mag 


0.039 mag 


Equation [3] 


0.15 dex 


0.035 dex 


0.021 dex 


0.017 dex 


Equation S] 


0.21 mag 


0.063 mag 


0.063 mag 


0.051 mag 


Equation [5] 


0.18 dex 


0.038 dex 


0.056 dex 


0.066 dex 


Equation |9] 


0.46 mag 


0.101 mag 


0.114 mag 


0.122 mag 


Equation [TOl 


0.39 dex 


0.092 dex 


0.063 dex 


0.105 dex 



Note. — The RMS difference A/ between our adopted form for Equations Ell5] and [MTUl and 
three alternate forms with comparable goodness-of-fit, averaged over the parameter space where 
each equation is valid. Were we to fit a different-degree polynomial for one of these equations, 
our predictions would be changed by roughly A/. Since A/ is typically much smaller than the 
residuals around our fits, the final choice of form for Equations [2]|5] and IMTOl does not affect our 
results. 
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Table 3. Mean Magnitude and Parallax in the Hipparcos Sample 
B-V logiJ^j^ < -5.0 -5.0<logi?^j^ <-4.5 -4.5 < logiijjj^ 



Mean V Magnitude 



0.5 


- 0.6 


6.45 ±0.14 


7.01 ± 0.04 


7.03 ± 0.09 


0.6 


- 0.7 


7.29 ±0.08 


7.60 ± 0.04 


7.67 ±0.08 


0.7 


- 0.8 


8.06 ± 0.09 


8.07 ±0.05 


8.32 ±0.09 


0.8 


- 0.9 


8.62 ± 0.11 


8.53 ±0.06 


8.42 ±0.12 


0.9 


- 1.0 


9.03 ±0.16 


9.03 ±0.06 


8.98 ±0.10 


Mean Parallax (mas) 


0.5 


- 0.6 


32.1 ±2.3 


29.1 ±0.7 


29.1 ± 1.2 


0.6 


- 0.7 


28.2 ± 1.0 


29.6 ±0.6 


30.2 ± 1.2 


0.7 


-0.8 


30.4 ± 1.4 


31.8 ± 1.1 


29.5 ± 1.5 


0.8 


- 0.9 


32.5 ± 1.6 


32.8 ± 1.3 


36.5 ±3.4 


0.9 


- 1.0 


30.2 ±2.0 


32.2 ±0.9 


32.1 ± 1.3 



Note. — Mean V magnitudes and parallaxes, with la errors on the moans, 
as a function of activity. All pairs of parallaxes are indistinguishable at 95% 
confidence. For B — V < 0.7, the least active bin is significantly brighter. 
However, since we see the strongest UV-activity correlations among the red- 
dest stars, the correlation of V with activity among the blue stars cannot be 
responsible for our trends between UV flux and activity. 
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Table 4. Predicted Activity Levels for Selected non-Hipparcos Stars 





Name 


FUV-V 
(mag) 


<^FUV 
(mag) 


NUV-V 
(mag) 


0"JV(7V 
(mag) 


B-V 
(mag) 


95% CI 
logR^K 


Observed 
logH^K 


"'°^«HK 




BD+01° 0306 






7.65 


0.01 


0.96 


(-5.10, -4.52) 


-5.07 


0.005 


CI* 


NGC 2632 KW 127 


11.17 


0.13 






0.60 


(-5.01, -4.31) 


-4.59 




CI* 


NGC 2632 KW 217 


9.93 


0.04 






0.51 


(-4.90, -4.20) 


-4.50 






HD 1342 


10.42 


0.14 


4.37 


0.00 


0.57 


(-4.84, -4.14) 


-4.89 


0.042 




HD 103195 






7.46 


0.01 


0.96 


(-5.06, -4.48) 


-4.87 






HD 115405 






6.33 


0.01 


0.85 


(-4.80, -4.22) 


-3.82 






HD 131157 


13.09 


0.64 


6.07 


0.009 


0.66 


(-5.51, -4.81) 


-3.64 


0.124 




HD 205724 






6.70 


0.03 


0.84 


(-5.15, -4.57) 


-4.71 


0.042 



Note. — 95% confidence intervals for activity levels predicted from UV data. The sample consists of stars with B — V 

and Ca H measurement s from iDuncan et al. ( 1991) that are not in the Hipparcos sample we used to calibrate the UV- 

activity relation. Where iDuncan et al. n^9l|) quoted multiple measurements, we selected one at random to ensure we did 

not suppress extreme or outlying measurements by averaging. Since the scatter of the measurements (column o-j^^g j^i ) is 

HK 

small, the choice of epoch does not affect our results. Most predictions are within the 0.15-0.18 dex uncertainty of Equations[3] 
andO but there are some exceptions, notably HD 131157. Aside from the two NGC 2632 members, none of these stars are 
discussed in the literature, so we have no information on whether they are interacting binaries, extreme metallicity stars, or 
otherwise unusual. 



